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AN ELASTIC ANALYSIS OF STRESSES IN 
A UNI AXIALLY LOADED SHEET CONTAINING AN 


INTERFERENCE -FIT BOLT 

By John H. Crews, Jr. 

Langley Research Center 

SUMMARY 

The stresses in a sheet with an interference -fit bolt have been calculated for two 
sheet -bolt interface conditions: a frictionless interface and a fixed (no-slip) interface. 
The stress distributions were calculated for various combinations of sheet and bolt 
moduli. The results show that for repeated loading the local stress range is significantly 
smaller if an interference bolt is used instead of a loosely fitting one. This reduction in 
local stress range is more pronounced when the ratio of bolt modulus to sheet modulus 
is large. The analysis also indicates that currently used standard values of interference 
cause yielding in the sheet. 


INTRODUCTION 

Laboratory studies have shown that structural joints with interference -fit bolts have 
longer fatigue lives than those with loosely fitted bolts (ref. 1). As a result, interference- 
fit bolts are used in aircraft structures (refs. 2 and 3), although the underlying phenom- 
enon which lengthens life is not completely understood. An analysis of the phenomenon 
was presented in reference 4, but unfortunately, the solution does not provide a clear 
explanation of the phenomenon nor does it lend itself to a systematic analysis of the prob- 
lem. The purpose of the present report is to present a more complete elastic analysis 
of stresses in a sheet containing an interference-fit bolt and to discuss the influence of 
interference fits on fatigue life. 

The analysis is based on a general solution presented in references 5 and 6 for a 
uniaxially loaded infinite sheet with a circular hole containing an elastic disk insert. Two 
special solutions were examined - the first was for a frictionless interface between the * 
sheet and bolt; the second was for a no-slip interface. These two idealizations bracket 
the real interface behavior, which is expected to involve some slip. 

Results are presented in terms of applied stress, interference, and material con- 
stants for the sheet and bolt. Stress distributions are presented for typical combinations 



of sheet and bolt materials and are used to explain the effects of interference bolts on 
fatigue life. Applicable ranges for the solution are presented in terms of the combina- 
tions of applied stress and interference that produce either yielding of the sheet or sepa- 
ration at the sheet-bolt interface. 

SYMBOLS 

The units used for the physical quantities defined in this paper are given in the 
International System of Units (SI) (ref. 7) and in the U.S. Customary Units. The measure- 
ments and calculations were made in the U.S. Customary Units. 

a^ag^, unknown coefficients in complex stress functions for a sheet with a circular 
hole 

b^bgjbg unknown coefficients in complex stress functions for a disk 
D diameter of hole, mm (in.) 

E Young's modulus, MN/m^ (ksi) 

I interference, or difference between bolt and hole diameters, mm (in.) 

R radial distance to bolt-sheet interface, mm (in.) 

r ,6 polar coordinates, mm (in.) and deg 

S nominal stress, MN/m^ (ksi) 

v r ,vg| displacements, mm (in.) 

x,y Cartesian coordinates, mm (in.) 

z complex variable, x + iy 

v Poisson's ratio 

" o rr ,<jQQ normal stress components, MN/m^ (ksi) 

o vv ,oqq normal stress components caused by interference fit between bolt and sheet, 
MN/m 2 (ksi) 
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yield stress, MN/m^ (ksi) 


r r g shear stress component, MN/m^ (ksi) 

0, 0 complex stress functions 


Subscripts: 

1 denotes sheet 

2 denotes bolt 

s separation 


ANALYSIS 

Figure 1 shows the problem analyzed in this study. The sheet was assumed to be 
infinite and the interference-fit bolt was represented by a disk equal in thickness to that 
of the sheet. The analysis was further simplified by assuming a state of plane elastic 
stress for both the sheet and bolt. This analysis becomes invalid when the hole deforms 
so much that the sheet separates from the bolt. 

The solution for this problem was taken from references 5 and 6 and is briefly 
described in the appendix. The stress equations for the sheet are as follows (see 
eqs. (A3) to (A5)): 



The undetermined coefficients a^ were evaluated for both the frictionless and no-slip 
cases by specifying stress and displacement conditions at the sheet -bolt interface, as 
explained in the appendix. 
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Because the analysis was based on the assumption of elastic material response, 
local yielding of the sheet material represents a limiting condition for the solutions. 
Yielding can be caused by high interference stresses alone or by the combined effect of 
interference stresses and stresses produced by remote loading of the sheet. Yielding 
was assumed to be governed by the Mises yield condition: 



°rr 1 (y ee 1 


-,2 

7 ee 1 


3t: 


rfli 


1/2 


( 2 ) 


Yielding is usually expected to initiate at the points of maximum stress concentration 
which lie on the interface boundary at the x-axis; however, for the no-slip case the inter- 
face shear stresses reach a maximum away from this axis and therefore can change the 
location for initial yielding. This possibility was examined by evaluating equation (2) at 
5° intervals along the interface. The desired limits for the solution were expressed in 
terms of interference and applied stress levels which cause yielding for various ratios 
of Young's moduli Eg^E^. 

In addition to the assumption of elastic material behavior, the sheet was assumed to 
be in contact with the bolt along the entire interface. Loading of the sheet eventually 
causes separation at the intersection of the y-axis and the interface. Since separation 
corresponds to zero-load transfer at these two points, a condition for the onset of sepa- 
ration was established by setting a rr ^R, 0 equal to zero. The result was an equation 

for the stress S s , which when applied to the sheet corresponds to incipient separation at 
the interface. 


For the frictionless case, the separation stress was found from equation (A18) to be 


S« = 


IEi 

2R 


Eo 


3 + "2 + ( 5 “ 


Ec 


9 - 5, 2 + (11 + 5^ 


(3) 


where I is interference, E is Young's modulus, v is Poisson's ratio, and the sub- 
scripts 1 and 2 denote the sheet and bolt, respectively. For the no-slip case, equa- 
tion (A25) gave 


Ss — 


IE 1 

~2R 


E, 


1 + v 2 + (3 - 


3 - + (5 + 


E r 

t 

Ei 


(4) 
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RESULTS AND DISCUSSION 


Because fatigue failures usually occur in the sheet rather than the bolt, only stresses 
in the sheet will be discussed. Stresses for a sheet with an interference-fit bolt are com- 
pared with those for a sheet with an open hole. For comparison, stresses for the open- 
hole case are presented in figure 2, where the local stresses have been normalized by 
the applied stress S. As expected, the normalized stress peaked at <*qq = 3 at 0=0. 
The familiar distributions of Oqq and a rr along the transverse axis of symmetry are 
shown in figure 2(b). Since only stresses for the sheet are presented, the subscript 1 
has been dropped. 


Frictionless Interface 

Stresses in figure 3 are for sheet and bolt materials with identical moduli. For con- 
venience, the stresses are normalized by S s , the applied stress for incipient separation 
at the interface. The interference stresses, produced solely by the interference fit, are 
indicated by a rr and Oqq and are represented by dashed lines; these stresses are 
uniform along the boundary of the hole (fig. 3(a)) and, as expected, approach zero away 
from the hole (fig. 3(b)). The combined effect of interference and the applied stress is 
shown by the solid curves. If the applied stress is cycled between 0 and S s , the stresses 
in the sheet cycle between the dashed and solid curves. Comparison of these curves in 
figure 3(a) shows that the maximum local stress range is only 1.75 (1.25 to 3.0), which is 
considerably smaller than the range 3.0 for the case of a loosely fitting bolt. This 
smaller local stress range is responsible for the improvement in fatigue life attributable 
to interference -fit fasteners. It should be noted that this smaller local stress range is 
accompanied by a larger local mean stress which tends to reduce the benefit of the inter- 
ference fit on fatigue life. 

E 2 

The stresses in figure 4 for — = 3 represent elastic behavior typical of an alumi- 

E 1 

num sheet with a steel bolt. The maximum value of o^q again equals 3.0, but since 
Oqq = 1.67 for this case, the maximum local stress range is only 1.33, compared with 
E o 

1.75 for — = 1, and 3.0 for the loose bolt. Thus, the local stress range is smaller for 

larger values of E 2 /Ej. This effect is illustrated in figure 5. As e 2 /e^ increases 
from zero, the local stress range decreases rapidly and approaches the limiting value of 
e 2 

about 0.95 for — = 00 (rigid bolt). For cyclic loading from 0 to S s , cycles 

E 1 y 
between values given by the dashed and solid curves. 

Limits for the frictionless-fit solution are presented in figure 6 in terms of curves 
corresponding to separation of the sheet -bolt interface and yielding of the sheet. The 
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separation curves, calculated from equation (3), show only a minor dependence on E 
In contrast, the yield curves, equation (2), are strongly dependent on For a 

given E 2 /E 2 value, the range of validity for the present solution is represented by the 
region under a corresponding pair of separation and yielding curves. In all frictionless- 
fit cases, yielding started at the intersection of the x-axis and the interface. 

Note that the intersection of the yield curves and the I/D axis represent yielding 

Eo 

due to interference alone. For the typical values of — = 1 and 3, yielding occurs for 

E 1 

0.0081 and 0.0062, respectively. Since standard values of I/D for interference 

bolts currently in use range from 0.007 to 0.017 (for D = 6.35 mm (0.25 in.), ref. 8), 
yielding due to interference -bolt installation is common. The less widely used "low 
interference" bolts (ref. 8) have I/D values of 0.0008 to 0.0060, which lie within the 
range of the present elastic solution. 


No-Slip Interface 

E 2 

Stresses in a sheet with a no-slip interference bolt are shown in figure 7 for — = 1. 

E 1 

The interference stresses for this case appear smaller than those for the corresponding 
frictionless case = 1.0 compared with 1.25 from fig. 3) because the value of S s 
used to normalize is larger for the no-slip case. This larger value of S g can be 

attributed to the shear stresses that develop at the no-slip interface when the sheet is 
loaded (fig. 7(a)). These stresses cause larger load transfer across the interface and 
therefore less load transfer around the bolt hole. Consequently, the strains (displace- 
ments) in the sheet that contribute to separation are smaller and thus permit larger 
values of S s . 


When stress is applied to the sheet, the maximum value of Oqq increases from 1.0 

Eo 

to 2.0 - a local stress range of 1.0. Since — = 1, and both shear and normal stresses 


are transferred across the interface, the hole is in effect eliminated by the presence of 
the bolt. Therefore, a local stress range equal to the applied stress range is expected. 
This is illustrated in figure 7(b), where the Oqq curve is identical with the curve 

when it is displaced vertically by 1.0. 

E 2 

Stresses for — = 3 are presented in figure 8. Although the interference stress at 
El / e 2 \ 

the interface is larger than that for the previous case — =11, the maximum local stress 


increases by only 0.45 to a value of 1.73 when the sheet is loaded. Thus, for — = 3, a 

El 1 
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no-slip interference fit is very effective in reducing the local cyclic stress range (0.45 
compared with 3.0 for the loosely fitting bolt). In fact, the local stress range at the bolt 
hole is smaller than that away from the hole, as shown in figure 8(b). 

The general influence of on local stress range for the no-slip case is shown 

in figure 9. Local stress range decreases abruptly as Ep/E t increases from zero, 

/Eo \ 

and it approaches zero for the rigid-bolt case — = °° . Note that for the loading used 

\ 1 / 

(0 to S s ), the local mean stress is independent of Eg/Ej. Therefore, the decrease in 
local stress range is achieved without the accompanying increase in local mean stress 
observed for the frictionless-fit case. 


Limits for the no-slip interference fit solution are shown in figure 10. These are 
similar to the curves for the frictionless fit in figure 6, but comparison shows that the 
range of applicability for the no-slip case is somewhat larger; that is, higher stresses 
can be applied without causing either separation or yielding. The point of initial yielding 
at the interface is shifted away from the x-axis by the shear stresses that develop along 

Eo o 

the interface when the sheet is loaded. For — = 1 and 3, this shift is less than 5 ; but 
„ E 1 
Eo o 

for — = °°, initial yielding occurs at about 9 = 25 . 

E 1 


CONCLUDING REMARKS 


The stresses in a uniaxially loaded sheet containing an interference-fit bolt have 
been calculated. These calculations explain why interference bolts lengthen fatigue life of 
aircraft structural components. Two special cases were examined - the first was based 
on a frictionless interface between the sheet and bolt; the second was based on a no-slip 
interface. These two idealizations bracket the real behavior at the interface, which is 
expected to involve some slip between the sheet and bolt. 

Both solutions show that for cyclic loading the range of local stresses at the bolt hole 
can be reduced significantly by using an interference -fit bolt instead of a loosely fitting 
bolt. The reduction is larger for the no-slip case and is more pronounced when a high- 
modulus bolt is used in a low-modulus sheet. 

For the frictionless case, the local stresses caused by interference and those caused 
by an applied stress S s , corresponding to incipient separation at the interface, combine 
to produce a maximum circumferential stress at the bolt hole equal to three times the 
applied stress. This is identical with the stress at a hole with a loosely fitting bolt. 
Therefore, the reduction in local stress range is due solely to the largely minimum local 
stress caused by the interference fit. 


7 



For the no- slip case, the maximum local stress is smaller than that for an open 
hole. The reduction in local stress range results from a smaller maximum local stress 
combined with a larger minimum stress due to interference. For cyclic loading from 0 
to S s , the local mean stress for this case is equal to that for the case of a loosely fitting 
bolt. 

An analysis of the stresses due solely to interference reveals that the standard 
values of interference currently in use usually cause sheet yielding during bolt installa- 
tion. Therefore, elastoplastic methods are required for a complete analysis of struc- 
tural members containing interference -fit bolts. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 12, 1972. 
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APPENDIX 


STRESS AND DISPLACEMENT EQUATIONS 

The desired stress and displacement equations for the problem of an interference- 
fit bolt in an infinite sheet under uniaxial load can be obtained from either reference 5 
or 6 . In reference 5, Muskhelishvili presented a solution for an unloaded infinite sheet 
with an interference-fit disk and a solution for a uniaxially loaded sheet with a noninter- 
ference disk as an insert. If the disk is assumed to represent a bolt, the desired equa- 
tions may be obtained by superposition of these two solutions. In reference 6 , Savin pre- 
sented a solution for a uniaxially loaded sheet with a hole reinforced by rings. He 
specialized this solution for the case of a single ring, then introduced interference, and 
finally set the inner radius of the ring equal to zero to represent an interference -fit disk. 
For this specialized case, the Savin solution reduces to that obtainable by superimposing 
the two Muskhelishvili solutions. 

The following development, although mathematically equivalent to the Muskhelishvili 
and Savin solutions, is presented in a more direct and complete manner. 

General Stress and Displacement Equations 

The stress and displacement equations were developed by the complex-variable 
approach (see, e.g., ref. 5). Accordingly, stresses a rr , o Qd , and T r0 and the dis- 
placements v r and v 0 were expressed in terms of complex stress functions (p( z) 
and i//(z). For the case of plane stress, 

°rr + a ee = 2 |^'( z ) + ^'( z )] 

°00 “ CT rr + 2iT r0 = 2 |j^"( z ) + ^'(zSJe 210 

(v r + iv 0 ) = <Mz) - z 0 '(z) - j//(z)|e " 10 

1+iA 1 + v J 

The stress functions c^(z) and ^(z) for the sheet were taken from reference 5. 
After a slight modification to change the loading from the x-direction to the y -direction, 
as shown in figure 1, they were expressed as 

SR, / „ R,\ 

0 X ( Z ) = _i ( — + a 2 — j (A2a) 
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APPENDIX - Continued 


SRJ „ R 1 R/ 

^»=-— -2_ + a 2 - + a 3 — 


(A2b) 


where R^ is the radius of the hole and a^, a 2 , and a^ are unknown coefficients. The 
functions <^(z) and j/^(z) were substituted in equations (Al) which were then solved 
for stresses and displacements; in terms of polar coordinates these equations were 
expressed as 


a 2/ R l 


-|l + 2ai^) +| a 3 (^) |cos20 


o / 2.9 /Ri \ q /R-i \ 

a 00 (r,0) = |a + — ( — j + 1 + 2 3\r ) C0S2e 


T r« 1 < r ' ,, > = ! 1 - a l(v-) ‘ f ^(t') 81,129 

v a-,) SR l( 1 ^l)/' 2 ( 1 - 1 J l) r , a 2^1 , J~ 2r , 4 a l R l 

r l ’ \ 1+ ^ Rj r r^l + 


+ a 3\“/ COS20 


v fl (r,0) = 




2r ^ “ v \ -Ri (RA 

— - 2a 1 — — — + aql — — ) sin 20 

Rj 1 1 + v 1 r r / 


The stress functions for the bolt were (ref. 5) 


. , SR 2 , / z \ 3 b 2 z 

2 4 ‘taj "2 


!^ 2 (Z) = - 


SR 2 /b 


4 VR2 
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APPENDIX - Continued 


where R 2 is the bolt radius and b^ are unknown coefficients. Substitution of 0 2 (z) 
and ^ ( z ) into equations (Al) yields 


W r ’ e> ° I 


bg 

+ — - cos 28 


(A9) 


and 



v fl (r , 8 ) 
°2 



(1 + r 2 )b 3 



sin 29 


(A10) 
(All) 
(A 12) 

(A 13) 


The general stress and displacement equations (A3) to (A7) and (A9) to (A13) were 
used in each of the two following solutions. 


Frictionless Interface 


An interference fit between a disk of radius R 2 
expressed as 


and a hole of radius R^ can be 


«2 - R l = |- >0 


where I is the interference. After the disk has been inserted into the hole, the disk- 
sheet interface lies at radius R, where R 1 < R < Rg. For simplicity, Rj and R 2 
are replaced by R in the stress and displacement equations; since I is much smaller 
than R^ and Rg, only small errors (less than 1 percent) are introduced. 

For the case of an interference fit with a frictionless interface, the following condi- 
tions must be satisfied on the interface: 

a rr (R,0)= a rr CM) (A 14) 

X 2 
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APPENDIX - Continued 


T V Q^ 6 ) = 0 

(A 15) 

T r 02 ( R ,0) = 0 

(A16) 

v (R, B) - v (R, 0) = \ 
r l r 2 2 

(A 17) 


The appropriate stress and displacement expressions from the general solution were sub- 
stituted into these interface equations. Two additional equations were obtained by sepa- 
rating algebraic and trigonometric terms in equations (A14) and (A16). The resulting 
set of six equations was solved for the six unknown coefficients in the general solution. 
The following coefficients were determined: 
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a ! = 


-2 + 


E, 

e. 

Ei 


Er 




4 1- 


ie i\e 2 


a 2 “ 2 - 


2SR/E 


1 


E, 


1 ' v 2 + (i + "l)^ 


Er 

L 

E, 


a 3 = 2 - 


3 + i>2 + (5 - 


Er 

L 

Ei 


-4 


Er 

£. 

E, 


b l = 


3 + 1^2 + (5 - 


E, 


b 2 = 


2 1 - 


ie i\e 2 


2SR/E 


1 " v 2 + ( X + "l) 


Er 

l 

E, 
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These coefficients were substituted into equations (A3) to (A5) to find the following equa- 
tions for stresses in the sheet: 
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APPENDIX - Continued 


and 


E 2 
6 =A 


^=|< 1 + 21 - 


E 


1 


3 + 1^2 + 


(5 - ,j) 


E 1 


4 * 1 

1- 1 

A 

d) 4 > 

E, 

3 + v 2 + ( 5 - 

j 


sin 26 


(A20) 


No-Slip Interface 

For the case of a no-slip interface the following interface conditions were used: 


rr 


r r 6 



(H,0) = a (R,0) 
1 rr 2 

(A21) 

(R,0) = T g (R,0) 
1 r0 2 

(A22) 

(R,6) - v (R, 6) = \ 
r 2 2 

(A23) 

(R,0) = v 0 (R, 9) 

(A24) 


Two additional equations were obtained by separating algebraic and trigonometric terms 
in equations (A22) and (A24). The following coefficients were determined from the set of 
six equations: 


Er 


E, 


a l “ 


-2 + 


1+v 2 + 


(3 - vj) 


E, 

i 

if 
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APPENDIX - Continued 


4 1 - 


a 2 - 2 


2SR/E! 


1 - ^2 + (! + ^i) 


E, 

i 

IT 


E, 

t 

IT 


a 3 - 2 - 


E, 


1 + v 2 + ( 3 - •’1)1: 


bj = 0 



2 



2SR/Ej 


E, 


1 - u 2 * ( l + *'i)e7 


and 


b 3 = 


E 2 

-8 

E 1 


1 + ^2 + ( 3 - ^i)i: 


These coefficients were substituted into equ a tions (A3) to (A5) to obtain the desired equa.- 
tions for stress in the sheet: 



(A25) 
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APPENDIX - Concluded 



(A26) 


(A27) 
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Figure 2.- Elastic stresses in a uniaxially loaded sheet with an open hole. 















(a) Stress distribution along the hole boundary. (b) Stress distribution on the transverse axis. 

Figure 7.- Elastic stresses in a loaded sheet with a no-slip interference -fit bolt. 
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— = 1; incipient separation. 







Incipient separation. 
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